Role of osmotic and hydrostatic pressures in bacteriophage genome ejection 
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A critical step in the bacteriophage hfe cycle is genome ejection into host bacteria. The ejection 
process for double-stranded DNA phages has been studied thoroughly in vitro, where after triggering 
with the cellular receptor the genome ejects into a buffer. The experimental data have been inter- 
preted in terms of the decrease in free energy of the densely packed DNA associated with genome 
ejection. Here we detail a simple model of genome ejection in terms of the hydrostatic and osmotic 
pressures inside the phage, a bacterium, and a buffer solution/culture medium. We argue that the 
hydrodynamic flow associated with the water movement from the buffer solution into the phage 
capsid and further drainage into the bacterial cytoplasm, driven by the osmotic gradient between 
the bacterial cytoplasm and culture medium, provides an alternative mechanism for phage genome 
ejection in vivo; the mechanism is perfectly consistent with phage genome ejection in vitro. 

PACS numbers: 87.15.A, 82.35.Pq, 82.39.Pj, 87.15.N 



I. INTRODUCTION 

Bacteriophages — phages for short — are viruses that 
infect bacteria. In order to infect a host, the phage tail 
(> 95% phage virions contain tails) attaches to a bac- 
terium, which triggers the opening of the tail, and the 
genome (usually double-stranded DNA) is delivered, in 
some cases together with protein molecules, into the cy- 
toplasm. Inside the capsid, the genome is packaged at 
near crystalline density (~500 mg/ml for many double- 
stranded DNA phages [1]), where it is in a state of high 
free energy due to DNA self-interactions (causing high 
osmotic pressure within a mature phage) and, to a lesser 
extent [5], the mechanical energy cost incurred in bend- 
ing the polymer into a compartment with a radius com- 
parable to DNA's persistence length. Energy for genome 
packaging is provided by the bacterial cell in which the 
phage was assembled. From a purely thermodynamic 
perspective, the stored free energy within the capsid can 
provide the driving force for DNA ejection, with the DNA 
acting essentially as a "loaded spring" that expands when 
opening of the tail is triggered. (There may be separate 
gates at the junction of the capsid and the tail tip, but 
we only consider the situation where all gates are either 
open or closed.) The loaded spring concept has been 
formalized in a continuum mechanics model |3H5] . Keep- 
ing in mind that the model applies only to equilibrium 
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states, and not to the dynamics or molecular mechanism 
of the process, it successfully explains in vitro data for 
phage A DNA ejection. Such experiments include the 
addition of capsid-permeable DNA-condensing agents or 
capsid-impermeable osmolytes to arrest ejection when 
the DNA remaining in the capsid reaches equilibrium 
with the buffer [BHS] and time-resolved light scattering 
[10] . In addition, fluorescence microscopy experiments 
following in vitro A DNA ejection in real time at the 
single-molecule level show that the genome is ejected in 
a single step at an ion-dependent rate, but one that is 
not a simple function of the amount of DNA remaining 
in the capsid [HI [H] . 

Experimental data on osmotic suppression of DNA 
ejection of phage SPPl [Hj are consistent with the con- 
tinuum mechanics model, but those on T5 are less so. 
The extent of T5 DNA ejected after tail triggering is not 
proportional to low external osmotic pressure and dif- 
ferent populations co-exist at equilibrium [13]. Further, 
DNA remaining in individual capsids was shown to un- 
dergo several phase transitions as ejection proceeds |14) . 
and there are distinct pauses where DNA exit from in- 
dividual particles is temporarily halted [TSl US . These 
pauses do not correspond to the natural nicks in the T5 
genome, which are also not required for pausing, and each 
re-starting event requires overcoming the same activation 
energy [TMT5] . It is not known what causes pausing dur- 
ing T5 DNA ejection. 

It is natural to imagine that the loaded spring concept 
extrapolates directly to DNA ejection in vivo, where the 
ejection medium consists of the bacterial cytoplasm in- 
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model that retains its key physical ingredients. Partic- 
ular emphasis is given to the fact that the hydrostatic 
pressure inside the viral capsid is free to adjust in re- 
sponse to DNA ejection, a factor that has been largely 
ignored. We show that such an analysis can assist the in- 
terpretation of existing experiments, as well as guide the 
creation of artificial model systems in which the physics 
of genome ejection can be further elucidated. 



FIG. 1: (a.) In vitro configuration involves capsid and envi- 
ronment compartments, each with its own osmotic (tt) and 
hydrostatic (P) pressure. As DNA is ejected (red arrow), its 
volume is replaced by water permeating through the capsid 
and flowing in through the tail (blue arrows), (b) In vivo 
conflguration with a third compartment, the interior of the 
bacterium. As the bacterium grows, it has a high internal 
osmotic pressure, vrt,, and pulls in water. The resulting water 
flow through the tail can have the opposite direction as in (a) 
and can drive DNA ejection, independently of the osmotic 
pressure inside the capsid. 

stead of a buffer solution. It has, however, recently been 
argued that many features of phage DNA ejection in vivo 
are incompatible with this simple concept [I9H2T] . These 
include the problems that the osmotic pressure inside 
growing bacteria should prevent complete ejection of the 
genome, the efficiency of phage infection being reduced 
when bacterial turgor is low (the stationary phase of bac- 
terial growth), ejection of proteins before (and also likely 
after) DNA ejection, and that some phages do not con- 
tain DNA packaged to high density. Furthermore, several 
phage types are known to eject their DNA slowly into 
the bacterial cytoplasm in two or more distinct phases. 
These observations suggest mechanism(s) other than the 
loaded spring may contribute to DNA ejection in vivo. 

Ejection into a bacterium differs in subtle ways from 
ejection into a buffer. In particular, ejection in vivo in- 
volves three distinct compartments (bacterium, capsid, 
environment) instead of two in vitro (capsid, environ- 
ment) . This is sketched in Fig. [ij Each compartment has 
distinct osmotic and hydrostatic pressures, with ramifi- 
cations for the transport of the genome out of the phage. 

Growing bacteria continuously draw in water from the 
environment. For this to occur, a pressure imbalance 
(osmotic and/or hydrostatic) must exist between the bac- 
terium and the environment. This imbalance represents 
a separate driving force for genome ejection, independent 
of the physics of the loaded spring model. The hydrody- 
namic model of genome ejection p!M2T] is consistent with 
many of the in vivo observations that are not readily ex- 
plained by the loaded spring concept [21]. We further 
argue below that the model has ramifications for the in- 
terpretation of in vitro data, in particular the dynamics 
of ejection. 

To date, the hydrodynamic model has remained rel- 
atively qualitative [201 HI] ■ The purpose of the present 
article is to present a simple explicit formulation of the 



II. DESCRIPTION OF THE MODEL 

In general, we consider three compartments, namely, 
the inside of the capsid, the bacterial cytoplasm and 
the surrounding culture medium (environment). We de- 
note the osmotic and the hydrostatic pressures of these 
compartments by {tTc, Pc}, {TTf,, Ph} and {tTe, Pe}, respec- 
tively. The variables n,, and Pe are treated as constants 
that characterize the environment; TTf, and Ph similarly 
define the state of the bacterial cytoplasm, tTc is a mono- 
tonically decreasing function of the amount of DNA re- 
maining in the capsid, and Pc is free to adjust to the dy- 
namics of DNA ejection. The thermodynamic potential 
for water movement from compartment 2 to compart- 
ment 1 is simply given by 

(1)12 = frl - Pi) -fr2 - P2) = -'^21- (1) 
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Assuming for simplicity that both water and the cap- 
sid are incompressible, during the ejection process the 
continuity equation for water movement implies that the 
rate at which DNA volume exits the capsid equals the 
volume rate W{t) of the net water entry into the cap- 
sid at all times. Water entry into (or exit out of) the 
capsid can occur via two distinct paths, namely, by per- 
meating through the capsid or by flowing through the 
tail tube. Having noted that the available cross-section 
of the tail-tube is much smaller than the capsid cross- 
section, we assume that permeation through the capsid 
is the dominant contribution to W{t). This simplification 
can easily be relaxed, but only at the cost of complicat- 
ing the analytical form of the final results without yield- 
ing additional insights into the underlying physics. We 
further assume that (i) ejection is sufficiently slow that 
all compartments (bacterium, capsid and environment) 
remains well-stirred (or equivalently, the water is fully 
equilibrated) in all compartments at all stages of ejection, 
i.e., the osmotic and hydrostatic pressures are uniquely, 
and uniformly, defined in all compartments, and (ii) wa- 
ter entry into the capsid from the environment is in the 
linear response regime. Under these conditions. 

Wit) = AdNA^DNa(0 « (f>ce{t). (2) 

Here ^dna is the cross-sectional area of the DNA, 
VuNA{t) is the speed of DNA ejection, and (j^ is the cap- 
sid's permeability to water. 
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Having neglected the rate of water flow through the tail 
tube in comparison to the rate of water entry into the 
capsid from the environment through the capsid shell, 
the dynamics of genome ejection both in vitro and in 
vivo are qualitatively the same; namely that water osmo- 
sis into the capsid from the environment is accompanied 
by DNA ejection from the capsid. The physical force 
responsible for DNA ejection is the hydrostatic pressure 
gradient acting on the DNA across the tail tube, and the 
corresponding work done in ejecting the DNA is derived 
from the free energy made available as water crosses the 
capsid shell [from the environment (with osmotic pres- 
sure TTg) into the capsid (with osmotic pressure tTc, where 

TTc > TTe)]. 



namic drag exerted by the fluid on the outgoing DNA, 
Fh = 2tt R()Lri(dq / dr)r=B.o, is 



Fh 



ln{R/Ro) 



2Rl 



2TTLr] 




^^DNA- (3) 



Note that with being a combination of fixed parameters 
and the freely adjustable capsid pressure, the rhs of Eq. 
([3]) contains two unknowns: the DNA velocity, VdnAj 
and the hydrostatic pressure of the capsid, Pc- In order 
to uniquely determine these quantities, we treat the in 
vitro and in vivo cases separately. 




tail 
wall 



FIG. 2: The flow velocity pattern of water through the gap 
between the DNA and the inner wall of the tail-tube, shown 
via a longitudinal cross-section. 

Although we neglect the contribution of the hydrody- 
namic flow of water in the tail for the purpose of de- 
termining the total flux of water into the capsid, flow 
in the tail does influence conditions locally. Specifically, 
any flow of water in the tail, v^, exerts a force on DNA, 
also located in the tail, via hydrodynamic drag. In or- 
der to estimate this force, we represent the DNA as a 
cylinder of radius Rq (i.e., Adna = t^Rq), and let the 
axis of the moving DNA coincide with the axis of the tail 
tube at all times during ejection. We denote the inner 
radius of the tail tube by R. We further apply no-slip 
boundary conditions at the surface of the DNA and at 
the inner surface of the tail tube. Under this simplified 
geometry, the velocity profile for water in the tail tube, 
q(r), exhibits cylindrical symmetry and is only a function 
of r, the radial distance measured from the axis of the 
DNA, and q{r) can be easily determined by solving the 
Stokes equation, r]r~^dr{rdrq{r)) = 4>/L (with no-slip 
boundary conditions). Here 77 is the (dynamic) viscosity 
of water, L is the length of the tail-tube and (j) represents 
(j)ce or (pcb for in-vitro and in vivo, respectively. The 
resulting flow velocity pattern of water through the gap 
between the DNA and the tail-tube, shown schematically 
in Fig. [2j The magnitude of the corresponding hydrody- 



III. HYDRODYNAMIC MODEL FOR IN VITRO 
EJECTION 

In the in vitro experiments the phage is immersed in 
a buffer solution, creating a two-compartment configura- 
tion: one compartment is the inside of the capsid, while 
the other is the environment (around the capsid). Before 
the opening of the tail is triggered, there is no net move- 
ment of water between these two compartments, conse- 
quently 0ce = and f2Tl 
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Pe + (tTc - TTe). 



(4) 



This is analogous to the classic case of the van 't Hoff 
water columns. 

Following triggering of the tail, DNA begins to exit 
the capsid while water simultaneously enters the capsid 
from the environment. In order to obtain an explicit 
expression for Vdna(Oi it is necessary to specify the re- 
lationship between Vdna(0 and the current state of the 
system. In principle, this requires knowledge of the free 
energy stored in the configuration of the DNA confined 
inside the capsid. This free energy is a sum of the in- 
teraction energy of the capsid content (DNA, water, ions 
and possibly other protein molecules), and the mechani- 
cal energy stored in the DNA due to its packaging in tight 
circles. The "mechanical" energy is however an order of 
magnitude smaller than the contribution from interac- 
tions [2] ; in the spirit of constructing the simplest model 
that captures the role of hydrodynamics, we therefore 
neglect it here. With this simplification, the force acting 
on the partially ejected DNA that experiences a hydro- 
static pressure difference Pce{t) = Pc{t) — Pe across the 
tail-tube is given by ylDNA-Pce(0- 

Thus, having denoted the effective mass of the mobile 
DNA by TO (where it is reasonable to assume that to is 
of the order of the mass of the DNA cylinder within the 
tail tube) and combining Eqs. ([1]), ([2]) and (H, we have 



mVBNAit) = Pce{t)ABNA- Fh{t), 



I.e. 



^DNA 



"l+^DNA , "2 



^DNA 



VbNA(i),(5) 
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and correspondingly, 



IV. HYDRODYNAMIC MODEL FOR IN VIVO 
EJECTION 



Pc{t) - Pe + T^ce{t) ~ ^^DNA (t) AdNa/(Tc 



(6) 



Equation ([s]) is easily solved as an integral equation for 
VdnaC^); in which m//3ADNA sets the (transient) time- 
scale beyond which the inertial term mVDNA(i) term is 
not relevant. The transient time-scale for typical tail- 
tube diameters and lengths is estimated to be in the range 
of microseconds. Given that the time-scales for in vitro 
ejection is in the range of seconds [TT], we neglect the 
intcrtial term, which leads us to 



VbNAW = - 



and Pc{t)=Pe- 



1- 



^DNA 



This result Q points to two distinct possible regimes 
for the dynamics of DNA ejection, (a) For high cap- 
sid permeability to water [i.e., (3ac ^ Adna], the DNA 
mobility through the tail tube is the rate-limiting mech- 
anism, in which case the hydrostatic pressure gradient 
Pc{t) — Pe simply follows the osmotic gradient ndt) — TTg 
in a manner analogous to Eq. Q. (b) In the opposite 
limit, namely for low capsid conductivity, water entry 
into the capsid is the rate-limiting step for DNA ejec- 
tion. In this scenario, the hydrostatic pressure gradient, 
Pc{t) — Pe, decreases abruptly when the DNA starts to 
eject following the triggering of the tail. This is the sce- 
nario considered for in vitro DNA ejection in Ref. [21) . 

Which limit, or where between the limits one should 
observe, is not a priori obvious. Given that the choice 
depends on the relative ease of water and DNA transport 
between the capsid and the environment, the choice likely 
depends on the phage studied. Some phages are resistant 
to severe osmotic downshifts; for example, T7 infects nor- 
mally when a suspension in ~ 5M CsCl is added directly 
to cells. Conversely, the classic Kleinschmidt picture of 
T2 DNA emanating from a ruptured capsid is a result of 
rapid dilution of a phage suspension in 2M ammonium 
acetate [55]. The resistance to osmotic downshift should 
correlate well to the mobility of water through the capsid 
membrane. In that case we would expect T7 to conform 
more closely to limit (a) than T2, which would be closer 
to limit (b). 

The final important point to note is that when all mo- 
tions of the DNA and water cease, i.e., at the end of 
ejection (partial or complete as it may be), then Eq. Q 
indicates that both osmotic and hydrostatic pressure gra- 
dients between the environment and the capsid become 
zero. These situations correspond to the osmotic sup- 
pression experiments [TH [13] ; complete suppression 
of ejection yields an estimate of the osmotic pressures in 
a mature phage particle: 25 atm for A, ~ 47 atm for 
SPPl, and ^ 16 atm for the deletion mutant T5st(0). 



For in vivo situations, it is important to appreciate 
that water uptake by the cytoplasm from the culture 
medium is very tightly regulated in a growing bacterium, 
which maintains an active osmotic gradient between the 
environment and the cytoplasm. A corresponding hy- 
drostatic pressure gradient (turgor) helps the bacterium 
enlarge [23] • For instance, the E. coli cytoplasm has a 
positive osmotic pressure of several atm above the envi- 
ronment (under various growth conditions the pressure 
has been estimated to vary between 2 and 15 atm, with 
3.5-5 atm being commonly accepted values (231[2i])- This 
osmotic pressure gradient is counter-balanced by a hydro- 
static pressure differential (turgor) ^ that enables the 
cell to enlarge during growth [26]. Gram-positive cells, 
such as Bacillus subtilis, the host for phage SPPl, have 
much higher turgor: 19 atm [27j . When the cell mem- 
brane is breached by an infecting phage, a new conduit 
opens up for water to enter from the environment into 
the cell cytoplasm by passing through the phage cap- 
sid. For the following analysis, we treat the cytoplasm 
as a uniform, well-stirred medium for which the osmotic 
and the hydrostatic pressures TTf, and Pb remain constant 
during ejection. This is reasonable as bacteria respond 
to osmotic perturbations within milliseconds. 

This situation is analogous to the in vitro case, with 
two important modifications: First, the thermodynamic 
potential difference for water transport between the envi- 
ronment (culture medium) and the bacterial cytoplasm, 
(/>fce, which, in the light of the above paragraph, is a 
strictly positive quantity for a growing bacterium, is given 
by 



'Jbe 



he + ' 



(8) 



Secondly, the hydrodynamics in the tail tube is driven 
by the themodynamic potential difference (pcti— ~4>bc) 
for water transport from the capsid into the bacterial 
cytoplasm, i.e.. 



Fh{t) = ai(l)cb{t) + a2VuNA- 



(9) 



We then use the continuity of water flow across the capsid 
shell ^ to write 



mVBNAit) = Pcb{t)A 



DNA 



Fh{t), 



I.e. 



wVDNA(i) , "i-I-Adna , 

^ 7rcb(t) H ^ 06e 



A 



DNA 



A 



DNA 



ai+^DNA , Oi2 



A 



DNA 



yDNA(t).(10) 



Once again, we neglect the inertial term (i.e., the tran- 
sient time m/ /3Ad-^a), and obtain the expressions for 
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T^DNA and Pc as 



Pc{t) = Pe 



A 



DNA 



T^ceXt) n <P6e-(ll) 



Equation (11) shows that the rate of DNA ejection is 



enhanced, irrespective of the osmotic pressure difference 
between the capsid and the bacterium; the driving force 
behind this enhancement being the hydrodynamic drag 
of water on the DNA in the tail tube, while the source 
of free energy being the osmotic imbalance between the 
bacterium and its environment |19H21j . As pointed out in 
Ref. [H], (j)be aiding the DNA ejection in vivo is superbly 
consistent with infection phenomenology. No phage can 
complete genome ejection without the 0f,e term in Eq. 
(Ill; in particular, phages related to P2, perhaps even P2 



and P4 themselves, cannot even initiate DNA release into 
the cytoplasm with only <f>cb at their disposal (a mature 
P4 phage has a very low internal osmotic pressure |28|). 
Note also that upon setting 0f,e = 0, we recover the in 
vitro results as of course we should. 



phage, a bacterium, and the surrounding medium. This 
model shows explicitly that the osmotic balance between 
a bacterium and its environment provides a second source 
of free energy for driving genome ejection that is indepen- 
dent of the free energy stored in the configuration of the 
DNA in the capsid. The model also demonstrates that 
water transport into the capsid and hydrodynamic drag 
on the DNA in the tail pipe can profoundly influence the 
rate of genome ejection, both in vitro and in vivo. Fur- 
ther refinements to the hydrodynamic model in the sim- 
ple form described here can be envisioned, for example 
accounting for a finite capsid elasticity or for the mechan- 
ical energy stored in the DNA, and this may prove nec- 
essary for a quantitative comparison with experiments. 
However, such refinements are not expected to signifi- 
cantly change the physical picture described here. 
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